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Abstract 

We consider discretizations of the hyper-singular integral operator on closed surfaces and show 
that the inverses of the corresponding system matrices can be approximated by blockwise low-rank 
matrices at an exponential rate in the block rank. We cover in particular the data-space format of 
4/,-matrices. We show the approximability result for two types of discretizations. The first one is a 
saddle point formulation, which incorporates the constraint of vanishing mean of the solution. The 
second discretization is based on a stabilized hyper-singular operator, which leads to symmetric 
positive definite matrices. In this latter setting, we also show that the hierarchical Cholesky factor¬ 
ization can be approximated at an exponential rate in the block rank. 


1 Introduction 


Boundary element method (BEM) are obtained as the discretizations of boundary boundary integral 
equations. These arise, for example, when elliptic partial differential equations are reformulated as in¬ 
tegral equations on the boundary T := dfl of a domain Q c K' / . A particular strength of these methods 
is that they can deal with unbounded exterior domains. Reformulating an equation posed in a vol¬ 
ume as one on its boundary brings about a significant reduction in complexity. However, the boundary 
integral operators are fully occupied, and this has sparked the development of various matrix compres¬ 
sion techniques. One possibility, which we will not pursue here, are wavelet compression techniques, 
[ Rat98RatOl ; ISch98al lvPSS971 lTau03l ' TW03 1. where sparsity of the system matrices results from the 
choice of basis. In the present work, we will consider data-sparse matrix formats that are based on block- 
wise low-rank matrices. These formats can be traced back to multipole expansions, [ Rok851 IGR97 1. 
panel clustering, I NI 188 HN891lHS93llSau92i) . and were then further developed in the mosaic-skeleton 
method, | Tyr00| , the adaptive cross approximation (ACA) method, l lBchOOl . and the hybrid cross ap¬ 
proximation (HCA), HBG05II . A fairly general framework for these techniques is given by the H- 
matrices, introduced in HHac991IGH031lGra01HHac09ll and the T/ 2 -matrices, IIHKSOOIlBdrlOalfBdrl Obi . 
Both %- and "^-matrices come with an (approximate) arithmetic and thus provide the possibility of 
(approximately) inverting or factorizing a BEM matrix; also algebraic approaches to the design of pre¬ 
conditioners for boundary element discretizations, both for positive and negative order operators, are 
available with this framework. Empirically, it has already been observed in IGraOlt iBeb05bll that such 
an approach works well in practice. 


Mathematically, the fundamental question in connection with the "H-matrix arithmetic is whether the 
desired result, i.e., the inverse (or a factorization such as an LU- or Cholesky factorization), can be 
represented accurately in near optimal complexity in this format. This question is answered in the 
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affirmative in the present work for discretizations of the hyper-singular integral operator associated with 
the Laplace operator. In previous work, we showed similar existence results for FEM discretizations 
IIFMPI 3bH and the discretization of the single layer operator, tfFMP13a i. Compared to the symmetric 
positive definite case of the single layer operator studied in [ FMP13a1 . the hyper-singular operator on 
closed surfaces has a one-dimensional kernel and is naturally treated as a (simple) saddle point problem. 
We show in Theorem l2.6l (cf. also Remark [2~71) that the inverse of the discretization of this saddle point 
formulation can be approximated by blockwise low-rank matrices at an exponential rate in the block 
rank. A corresponding approximation result for the discretized version of the stabilized hyper-singular 
operator follows then fairly easily in Corollarv l5.ll The approximation result Theorem l2.6l also underlies 
our proof that the hierarchical Cholesky factorization of the stabilized hyper-singular operator admits 
an efficient representation in the H -matrix format (Theorem 16. II) . 

The approximability problem for the inverses of Galerkin BEM-matrices has previously only been 
studied in IIFMPI3al for the single layer operator. In a FEM context, works prior to lFMP13bl include 
I Bl 103 Beb05a . Beb05b l. I ISch06l l. and HBorlOat . These works differ from I FMP13bl EMP13a l and 
the present paper in an important technical aspect: while IIFMPI3bl IFMP13all and the present analy¬ 
sis analyze the discretized operators and show exponential convergence in the block rank, the above 
mentioned works study first low-rank approximations on the continuous level and transfer these to the 
discrete level in a final projection step. Therefore, they achieve exponential convergence in the block 
rank up to this projection error, which is related to the discretization error. 

The paper is structured as follows. In the interest of readability, we have collected the main result 
concerning the approximability of the inverse of the discretization of the saddle point formulation in 
Section [2] The mathematical core is found in Section 0 where we study how well solutions of the 
(discretized) hyper-singular integral equation can be approximated from low-dimensional spaces (The¬ 
orem [3j]). In contrast to IIFMPI3afl . which considered only lowest-order discretization, we consider 
here arbitrary fixed-order discretizations. The approximation result of Section[3]can be translated to the 
matrix level, which is done in Section 0] Section [5] shows how the results for the saddle point formu¬ 
lation imply corresponding ones for the stabilized hyper-singular operator. Finally, Section [6] provides 
the existence of an approximate %-Cholesky decomposition. We close with numerical examples in 
Section [7] 

We use standard integer order Sobolev spaces and the fractional order Sobolev spaces H l / 2 (T) and 
its dual iF~ 1 / 2 (r) as defined in, e.g., I ISSI 11 . The notation < abbreviates < up to a constant C > 0 
that depends only on the domain fi, the spatial dimension d, the polynomial degree p, and the 7- 
shape regularity of 7) t . It does not, however, depend on critical parameters such as the mesh size h, the 
dimension of the finite dimensional BEM space, or the block rank employed. Moreover, we use ~ to 
indicate that both estimates < and > hold. 


2 Main Result 

2.1 Notation and setting 

Throughout this paper, we assume that C K'b d G {2, 3} is a bounded Lipschitz domain such that 
T := dQ is polygonal (for d = 2) or polyhedral (for d = 3). We assume that F is connected. 

We consider the hyper-singular integral operator W € L(iF 1 / 2 (r), H~ 1 ^ 2 ( T) ) given by 

Wv{x) = - 77 LO Kv)(x) = - 7 M^( 7 i n yG(x - y))v{y)ds y , iGT, 
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where G(x) = — ipz log |xj for d = 2 and G(x) = for d = 3 is the fundamental solution 

associated with the Laplacian. Here, the double layer potential K £ L(H 1 / 2 ( r), is given by 

Kv{x) := f r (-ff y G(x — y))v{y)ds y , where 7 “^ denotes the interior conormal derivative at the point 
2 € T, i.e., with the normal vector n(z) at z € T pointing into Q c and some sufficiently smooth function 
u defined in Q one requires 7 f z u = Vu(z) ■ n(z). 

The hyper-singular integral operator W is symmetric, positive semidefinite on H l/ 2 (T). Since T is 
connected, W has a one-dimensional kernel given by the constant functions. In order to deal with this 
kernel, we can either use factor spaces, stabilize the operator, or study a saddle point formulation. In the 
following, we will employ the latter by adding the side constraint of vanishing mean. In Section [5] we 
will very briefly study the case of the stabilized operator, and our analysis of Cholesky factorizations in 
Section [ 6 ] will be performed for the stabilized operator. 

With the bilinear form b(v,p) := p J r vds x , we get the saddle point formulation of the boundary 
integral equation 

Wcp = f on T 

with arbitrary / £ P[~ l / 2 (T) as finding (cp, A) £ Tf 1 / 2 (T) x R such that 

(Wc/>^)+b^ 7 X) = (f,^) \hp £ H 1/2 (T), (2.1a) 

b{(j),n) = 0 V 7 £ R. (2.1b) 

By classical saddle-point theory, this problem has a unique solution (<p, A) £ // 1 /2 (I') x R, since the 
bilinear form b satisfies an inf-sup condition, and the bilinear form (Wcp, ip) is coercive on the kernel 
of b(-, A), which is just the one-dimensional space of constant functions (see, e.g., itSS 1 ill ). 

For the discretization, we assume that T is triangulated by a (globally) quasiuniform mesh Th = 
{Tj,..., Tm } of mesh width h := max^ g T h diam(Tj). The elements Tj £ T, are open line seg¬ 
ments (d = 2) or triangles (d = 3). Additionally, we assume that the mesh 7/, is regular in the sense 
of Ciarlet and 7 -shape regular in the sense that for d = 2 the quotient of the diameters of neighbor¬ 
ing elements is bounded by 7 and for d = 3 we have diam(Tj ) < 7 I ?)| 1/2 for all Tj £ Th, where 
T ;) | = area(T_,) denotes the length/area of the element Tj. 

We consider the Galerkin discretization of W by continuous, piecewise polynomial functions of fixed 
degree p > 1 in S p,1 (Th ) := {u € C(T) : u\t € P p {T)\/T £ Th}, where P p (T ) denotes the space of 
polynomials of maximal degree p on the triangle T. We choose a basis of S p ' 1 (T ), which is denoted 
by Bh '■= {ipj '■ j = 1, • • •, N}. Given that our results are formulated for matrices, assumptions on the 
basis Bh need to be imposed. For the isomorphism $ : R^ —> S p,1 (Th), x i-£ ^ 7=1 we require 

h {d - 1)/2 ||x || 2 < ||$(x)|| i2(r) < h^' 2 ||x || 2 Vx £ R^. (2.2) 


Remark 2.1 The standard basis for p = 1 consists of the classical hat functions satisfying ipj(xi) = Sij 
and for p > 2 we refer to, e.g., BSch98bl|KS99[[DKP + 08ll . These bases satisfy assumption (12.21) . ■ 


The discrete variational problem is given by finding {f>h, A/,) € S p ' 1 (T) x R such that 

(Wcp h ,ip h )+b(ip h ,X h ) = ( f,tp h ) Vip h € S P,1 (T), (2.3) 

b{(ph,p) = 0 V/x £ R. 


Since the bilinear form b trivially satisfies a discrete inf-sup condition, the discrete problem is uniquely 
solvable as well, and one has the stability bounds 

Il0ti||/r i /2(r) + |A| < G7 ||/||#-i/2 ( r ), (2.4) 
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for a constant C > 0 which depends only on T. For / € L 2 ( T) and the L 2 -projection II L “ : L 2 (T) —>• 
S p,1 (Th), one even has the following estimate 


ll^ll/t 1 /2(r) + l-^l < C 


n L / 


l 2 ( r) 


<C\\f\\ L 2 {r) 


With the basis Bh, the left-hand side of (12.31) leads to the invertible block matrix 

W: =fW B 

\B t 0 


(2.5) 


( 2 . 6 ) 


where the matrix W e 


pNxN 


and the vector B € 


pTVxl 


are given by 


w jk = {Wipk,i>j), B j = (i/jj, 1), e B h - (2-7) 


2.2 Approximation of W 1 by blockwise low-rank matrices 

Our goal is to approximate the inverse matrix W 1 by 77-matrices, which are based on the concept that 
certain ’admissible’ blocks can be approximated by low-rank factorizations. The following definition 
specifies for which blocks such a factorization can be derived. 

Definition 2.2 (bounding boxes and //-admissibility) A cluster r is a subset of the index set X = 
{1,... , N}. For a cluster r CX, we say that B Rt C M. d is a bounding box if: 

(i) B Rt is a hyper cube with side length R r , 

(ii) snpp fi C B Rr for all i £ r. 

For an admissibility parameter // > 0, a pair of clusters (r, o) with r, o C X is //-admissible if there 
exist bounding boxes Br t , Bp> a satisfying ( 0 -(O such that 

min{diam(B Rr ), diam(B R(T )} < r/ dist(B Rr , B Ra ). (2.8) 

Definition 2.3 (blockwise rank-r matrices) Let P be a partition ofX x X and // > 0. A matrix W ^ G 
1 R' A ' x A ' is said to be a blockwise rank-r matrix, if for every rj-admissible cluster pair (r, a) € P, the 
block W%| rX(T is a rank-r matrix, i.e., it has the form W%| rX(T = ~S. ra 'Y^ (J with X rcr € Ml T l xr and 
Y rcr € l^l^l xr '. Here and below, |cr| denotes the cardinality of a finite set o. 

Definition 2.4 (cluster tree) A cluster tree with leaf size m ea f € N is a binary tree Tj with root X such 
that for each cluster r € Tj the following dichotomy holds: either r is a leaf of the tree and |r| < ni ea f, 
or there exist sons t' , t" € Tj, which are disjoint subsets of r with r = t' U t". The level function 
level : Tj —y No is inductively defined by level(X) = 0 and level(r 7 ) := level(r) + 1 for t' a son ofr. 
The depth of a cluster tree is depth(Tj) := max^-g^x level(r). 

Definition 2.5 (far field, near field, and sparsity constant) A partition P ofX x X is said to be based 
on the cluster tree Tj, if P C Tj x Tj. For such a partition P and a fixed admissibility parameter 
>1 > 0 , we define the far field and the near field as 

Bf ar := {(r,cr) € P : (r,or) is q-admissible}, P near := P\ P far . (2.9) 

The sparsity constant C sp of such a partition was introduced in hGraOl 1/ as 

C sp := max < max |{cr £ Tx : r x a G Pf ar }|, max |{r G Tj : r x a € Pf ar }| > . (2.10) 

IrgTx ctGTx I 
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The following theorem is the main result of this paper. It states that the inverse matrix W 1 can be 
approximated by an %-matrix, where the approximation error in the spectral norm converges exponen¬ 
tially in the block rank. 


Theorem 2.6 Fix an admissibility parameter // > 0. Let a partition P ofT x X he based on the cluster 
tree Tj. Then, there exists a blockwise rank-r matrix Vp such that 

||W-Vxiv - V W || 2 < C' apx C' sp depth(T I )iV( 2d - 1 )/( M - 2 )e- 6 " 1/(d+1) . 

The constant C apx depends only on f l, d, p, and the y-shape regularity of the quasiuniform triangulation 
7/,, while the constant b > 0 additionally depends on q. 


Remark 2.7 (approximation of inverse of full system) The previous theorem provides an approxi¬ 
mation Vp to the first N x iV-subblock V of the matrix W 1 = (pT ^ j - Since P € M iVxl 


is a vector, the matrix \p 
satisfying 




is a blockwise rank-r approximation to the matrix W 1 


W ” 1 




< C apx C sp depth(T x )^( 2d - 1 )/( 2 ' i - 2 )e- brl/([i+1) . 


Remark 2.8 (relative errors) In order to derive a bound for the relative error, we need an estimate on 
||W|| 2 > since < ||W|| 2 - Since W is symmetric it suffices to estimate the Rayleigh quotient. 

The continuity of the hyper-singular integral operator as well as an inverse inequality, see Lemma 13.51 
below, and (12.21 ) imply 


W 


X ’ l A 


< 

r^j 


1^11 z/ 1 / 2 (r) ^ 1^ ( v i 1) I 


< 

r^j 


h 1 IMli 2 (r) + 1^1 ll^llz, 2 (r) ~ h d 2 

Using h — we get a bound for the relative error 


W 1 - \p 


W 


1 < C apx C sp iV( d+1 )/( M - 2 )depth(T x )e- fer ' 1/(li+1) . 


( 2 . 11 ) 


3 Approximation of the potential 

In order to approximate the inverse matrix W 1 by a blockwise low-rank matrix, we will analyze how 
well the solution of (12.31 ) can be approximated from low dimensional spaces. 
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Solving the problem (12.31) is equivalent to solving the linear system 


B w ^ ?) Q - (5) 

with W, B from (12.71) and b G IR: V defined by bj = (/, ipj). 

The solution vector x is linked to the Galerkin solution (fi, from (12.31) via (fi, = J2 f=i ^j^j- 

In this section, we will repeatedly use the L 2 (T)-orthogonal projection 11 1,2 : L 2 (Y) —> S p,1 (Th ) onto 
S p,1 (Th), which, we recall, is defined by 

(u La v^ h ) = (vM W h eS p '\T h ). (3.2) 

The following theorem is the main result of this section; it states that for an admissible block (r, a), 
there exists a low dimensional approximation space such that the restriction to B Rt D T of the Galerkin 
solution <ph can be approximated well from it as soon as the right-hand side / has support in B Ra n T. 

Theorem 3.1 Let (r, a) be a cluster pair with bounding boxes B Rt , B Ra (cf Definition 12,21) . Assume 
77 dist (Br t ,B R( j ) > d\&m(B Rr ) for some admissibility parameter rj > 0. Fix q G (0,1). Then, for 
each k G N there exists a space Wk C S p,1 (7h) with dim Wk < C uim (2 + r/) d q~ d k d+1 such that for 
arbitrary f G L 2 (T) with supp / C B Ra D T, the solution fh of (12.31) satisfies 

^min ||^ h - w\\ L 2 ( BRtnr) < C hox h~ 1/2 q k \\n L ~ f\\ L 2 {r) < C box h~ 1/2 q k \\f\\ L 2 (r) . (3.3) 

The constants C<}■„„, C\, ox > 0 depend only on LI, d, p, and the "/-shape regularity of the quasiuniform 
triangulation Th- 

The proof of Theorem 13. II will be given at the end of this section. Its main ingredients can be summa¬ 
rized as follows: First, the double-layer potential 

u(x) := Kfi h (x) = J 7 " %G(x - y)fi h (y)ds y , x G R d \ T, 

generated by the solution fh of (12.31) is harmonic on f l as well as on Ll c := W l \ LI and satisfies the 
jump conditions 

[7ow] := 7o Xt w - 7o m w = 4>h € H 1/2 (T ), 

[d n u] := 7i Xt w-7l nt w = 0GiT- 1 / 2 (r). (3.4) 

Here, 7 o xt , 7 o nt denote the exterior and interior trace operator and 7® xt ,7j nt the exterior and interior 
conormal derivative, see, e.g., I ISS I il l. Hence, the potential a is in a space of piecewise harmonic 
functions, where the jump across the boundary is a continuous piecewise polynomial of degree p, 
and the jump of the normal derivative vanishes. These properties will characterize the spaces l~Lh{D) 
to be introduced below. The second observation is an orthogonality condition on admissible blocks 
(r, < 7 ). For right-hand sides / with supp / C Bj ia n F, equation (12.31) . the admissibility condition, and 
W = —7 imply 

- + A (VTi, 1) = 0 \/fi h G S p,1 (Th) with supp^fe C B Rt D T. (3.5) 
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For a cluster p C X, we define r p C T as an open polygonal manifold given by 


:= interior 


|J SUpp fij 


U'GP 


(3-6) 


Let D be an open set and D 
harmonic, if 


:= D nil, D + := Dn Q c . A function v G H 1 (D\ T) is called piecewise 



Vn • \7tpdx 


0 V^eCo 00 ^). 


Definition 3.2 Let I) C M f/ be open. The restrictions of the interior and exterior trace operators y("‘, 
to DFT are operators 7 o" f |r>nr : H l (D~) —> Lf oc (DFT) andy e 0 xt \Dr\V ■ H 1 (D + ) -4 Lf oc (DnT) 
defined in the following way: For any (relative) compact JJ C I) FI T, one selects a cut-off function 
p G Cq°(D) with p = 1 on U. Since u G H 1 (D~) implies pu G H l {Q), we have Jq lt pu € H l / 2 (T) 
and thus its restriction to U is a well-defined function in I 2 (U). It is easy to see that the values on U 
do not depend on the choice ofp. The operator yjfi 1 onr is defined completely analogously. 

In order to define the restriction of the normal derivative of a piecewise harmonic function v G // 1 (D \ 
r), let p G C°°(M d ) with suppp C D and p = 1 on a compact set U C D. Then, the exterior normal 
derivative d n {pv ) is well defined as afunctional in H~ l / 2 ( r), and we define d n v\jj as the functional 

{d n v\u,ip) = (d n (pv),ip), \/<p£H 1/2 (T) , supp p C U. 

Again, this definition does not depend on the choice ofp as long as p = 1 on U. 


Definition 3.3 For a piecewise harmonic function v G II 1 ( I) \ F), we define the jump of the normal 
derivative [Sn'uJI.Dnr on D FT as the functional 

([dnt>]|Dnr,¥>) := f Vv-Vipdx Vip £ Hq(D). (3.7) 

Jd+ud- 

We note that the value ([3 n ^]|Dnr, <p) depends only on p\Dnv in the sense that ([<9 n 'u]|.Dnr, p) = 0 
for all tp G Cy(D) with ip\Dnv = 0. Moreover, if [9 n n]|i)nr i s a function in L 2 (D n T), then it is 
unique. The definition (13.71) is consistent with (13.41) in the following sense: For a potential Kfh with 
fh F S p,1 (Th), we have the jump condition [ d n K<fh]\Drr = 0. 

With these observations, we can define the space 

Flh(D) := {v £ H 1 (D\T): v is piecewise harmonic, [<9n^]|Dnr = 0, 

3v G S p,1 (Th ) s.t. [ 7 on]|r)nr = ^|z)nr}. 

The potential u = Kfh indeed satisfies u G 'H/fl)) for any domain I): we will later take 1) to be a 
box Bj j. 

For a box Br with side length R, we introduce the following norm on H 1 (Br \ F) 

:= ll Vu lli 2 (s fl \r) + ^2 IHli 2 (Bfl\r) > 

which is, for fixed h, equivalent to the H 1 (Br \ r)-norm. 
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A main tool in our proofs is the nodal interpolation operator //, : C(T) —>• S p,1 (Th)- Since p + 1 > 
(recall: d G {2,3}), the interpolation operator //, has the following local approximation property 
for continuous, Th -piecewise iF p+1 -functions u G C(T) n -ffpw 1 (7/ l ) := {u G C'(r) : u\t G 
H P+1 (T)VT G Th}- 

||tt - < C'/i 2(p+1_m) |rt|^- P +i( T ), 0 < m < p + 1. (3.8) 

The constant C > 0 depends only on 7 -shape regularity of the quasiuniform triangulation Th, the 
dimension d, and the polynomial degree p. 

In the following, we will approximate the Galerkin solution on certain nested boxes, which are concen¬ 
tric according to the following definition. 

Definition 3.4 Two (open) boxes Ii r, B /,>/ are said to be concentric boxes with side lengths R and R' , 
if they have the same barycenter and Br can be obtained by a stretching of Br/ by the factor R/R' 
taking their common barycenter as the origin. 

The following lemma states two classical inverse inequalities for functions in S p,1 (Th), which will 
repeatedly be used in this section. For a proof we refer to IIGHS051 Theorem 3.2] and [ SSI 11 Theorem 
4.4.2], 

Lemma 3.5 There is a constant C > 0 depending only on LI, d, p, and the 7 -shape regularity of the 
quasiuniform triangulation Th such that for all s G [0,1] the inverse inequality 

IMIff-s(r) < ChT 1 / 2 |M| L 2 (d Vu G S p,1 (Th) (3.9) 

holds. Furthermore, for 0 < rn < £ the inverse estimate 

\\v\\ H t {T) < Ch m ~ l \\v\\ Hrn(T) , Vv G V P (T) (3.10) 

holds for all T G Th, where the constant C > 0 depends only on (l.p, (l and the 7 -shape regularity of 
the quasiuniform triangulation Th- 

The following lemma shows that for piecewise harmonic functions, the restriction of the normal deriva¬ 
tive is a function in L 2 on a smaller box, and provides an estimate of the L 2 -norm of the normal 
derivative. 

Lemma 3.6 Let 6 G (0,1), R G (0, 2 diam(fi)) be such that ^ < |, and let /iGi. Let Br, B( 1+ $)r 
be two concentric boxes of side lengths R and (1 + 6)11. Then, there exists a constant C > 0 depending 
only on (>, d. p, and the y-sliape regularity of the quasiuniform triangulation Th, such that for all 
v G U h {B( 1+l 5)H ) we have 

ll^» u llL2(B fln r) ^ Ch V2 ^ll Vu llL 2 (B(i +5)i j\r) + ll v lli 2 (B ( i + 7fl\r)) ■ ( 3 - n ) 


Proof: 

1 . step: Let 77 G FF 1 ’ 0O (M a! ) satisfy 0 < 77 < 1 , 7 = 1 on -B(i+,5/2)j?, suppr/ C B^+^r, and 
||Vjyll^oo]^ a ) < jp- In order to shorten the proof, we assume 7 J ^77 = 7 o xt ?7 € S' 1 , 1 ( Th) so 
that inverse inequalities are applicable. We mention in passing that this simplification could be avoided 






by using “super-approximation”, a technique that goes back to BNS74I (cf., e.g., I Wah91l Assump¬ 
tion 7.1]). Let us briefly indicate, how the assumption rj € S ,1,1 ( Th) can be ensured: Start from a 
smooth cut-off function rj £ Cq c (M' / ) with the desired support properties. Then, the piecewise linear 
interpolant lj L ?j € 5' 1,1 ( Th) has the desired properties on T. It therefore suffices to construct a suitable 
lifting. This is achieved with the lifting operator described in llSte70l Chap. VI, Thm. 3] and afterwards 
a multiplication by a suitable cut-off function again. 

2. step: Let 2 := K('-/" j " r)[v] ). Then with the jump conditions 

[d n z] = 0, [ 2 ] = 7o nt? 7M 


and the fact that v is piecewise harmonic, we get that the function v—z is harmonic in the box B(i + $ / 2 )r- 
Thus, the function w := V(n — z) is harmonic in B^ 1+S / 2 )r as well. It therefore satisfies the interior 
regularity (Caccioppoli) estimate 

ll Vu ’llL 2 (B (1+ j /4 ) i? \r) < Jx IHl£ 2 (B ( i+*/2)fl\r) 5 (3 ‘ 12) 

a short proof of this Caccioppoli inequality can be found, for example, in | [BH03i l. 

We will need a second smooth cut-off function rj with 0 < rj < 1, rj = 1 on Br, and suppp C 
B(i+s/A)R and ||Vry|| l°°{B( 1+s/4)r ) < Jr ■ The multiplicative trace inequality, see, e.g., [ BS021 , implies 
together with (13.121) and SR < 2diam(Q) due to the assumptions on 5 , R that 

WvMl l^b R nr) ^ llwllL2(B C i +V 4) fl \r) + WvM\L^(B ( 1 +s/ 4 )R \r) W V< T w )WL^B ( 1 +s/ 4 )R \r) 

~ ll^llL 2 (s ( i + 5 / 4) fl\r) + \^ w hz(B (1+5/4)R \r) IML 2 (i? (1+ ,, /4)K \r) + H Vu; lli 2 (B(i +5 /4)n\r) 

< 3 11 11 2 

~ SR " wl ' L2 (B (1+s ,2)RXr) ■ 

Therefore and with d n z = d n K(T(j' y i'j[i:}) = —ll / ( 7 ( l | n, 7 /[■(■]), we can estimate the normal derivative of 
v by 


W^n^W l 2 (b R nr) — 11^ ' n llL 2 (s H nr) + H^ z llL 2 (s H nr) 

~ -jjjpj H u; ll L2 (' B(i +5 /2) fi\ r ) + II^C 70 7 ?H)IL 2(i?J?nr ). 


Since the hyper-singular integral operator is a continuous mapping from H 1 (I j to L 2 (V) and the double 
layer potential is continuous from H l / 2 (T) to (see, e.g., ISS111 Remark 3.1.18.]), we get with 

h < 5R, the inverse inequality (13.91) (note that ( 7 o nt 7 )[n] is a piecewise polynomial), and the trace 
inequality 


\\dn v \\L 2 (B R nr) 


~ ^ w ^^(B ( 1 +s/ 2 )R \r) + ||( 7 o 7 7 )II 

~ -jjfpj (ll V ' u ll L2 ( B ( i +5/2)j?\ r ) + H V "lli 2 (B(i+ 72 )fl\r)) + h Ml(To v)[v}\\ H i, 2(r) 

< h - 1/2 (11Vn|| L 2 ( s (1+4/2)H \ r) + ||(7o n S)HLi/ 2(r) ) 

< h 1/2 (11Vn11 L 2 (B (1+4/2) H\ r) + ll ? ? ?; ll// 1 (B (1+i)J? \r)) 

< h 1/2 ^||Vu|| i2(B(1+s)iJ \ r) + IHlL 2 (B {1+i)J? \r)) > 
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which finishes the proof. 


□ 


The previous lemma implies that for functions in 'Hh{B(i+5)R)i the normal derivative is a function in 
L 2 ( B[{ nT). Together with the orthogonality properties that we have identified in (13.51) . this is captured 
by the following affine space TLh 0 (D, T p , p): 

Tih,o(D,r p ,/j,) := R h (D)n{v e H L (D\T): supp[j 0 v}\ Dnr CT p , (3.13) 

(d n v\Drr,^h) ~ fJ'i'iph, 1) = OVVh, G S p ’ l (Th) with supp^h C finT p }. 

Lemma 3.7 The spaces TLh(D) and 'Hh,o(D, T p , p) are closed subspaces of H l (D \ T). 

Proof: Let (i4) je pj C 'Hh(D) be a sequence converging to v G H 1 ( I) \ T). With the definition of 
the jump [ 7 o'c J ]|Dnr and the continuity of the trace operator from // 1 (Q) to L 2 (T), we get that the 
sequence [ 7 oU' J ']|.Dnr converges in L 2 oc (D n T) to [ 7 ov]|.Dnr 5 and since S p,1 (Th ) is finite dimensional, 
we get that [ 70 ^]|z>nr = v\dc\T with a function v G S p,1 (Th)- 

Moreover, for ip G Cq 2 (D ± ) we have 

(Vv,Vip) L2{D \ r) =lunJVvfVp) L2{DV) = 0 , 

so v is piecewise harmonic on I)\Y. By definition (13.71) and the same argument, we get \d n v\ \ onv = 0, 
and therefore TLh{D) is closed. The space 'Hh,o{D, F p , p) is closed, since the intersection of closed 
spaces is closed. □ 

A key ingredient of the proof of Theorem 13. ll is a Caccioppoli-type interior regularity estimate, which 
is proved by use of the orthogonality property (13.51) . 

Lemma 3.8 Let 5 G (0,1), R G (0, 2diam(fi)) such that < | and let T p C T be of the form 
dm Let Br, be two concentric boxes and let p G R. Then, there exists a constant C > 0 

depending only on LI, d, p, and the y-shape regularity of the quasiuniform triangulation 7j, such that 
for all v G 'Hh,o{B(i +S )R,T p , p) 

ll Vt, llL 2 (s fl \r) - c + + ^) R )^ d 1)/2 • (3.14) 

Proof: Let q G iT 1 (M a! ) be a cut-off function with supp?] C Bn + g/ 2 )R, q = 1 on Br, and 
ll^ 7 ?llL°°(B (1+4 - )i j) ~ Jr- in the proof of Lemma lT6l we may additionally assume that = Jq xt q 
is a piecewise polynomial of degree 1 on each connected component of T D B(\ + s)R- Since h is the 
maximal element diameter, 8 h < SR implies T C Ifq+npi for all T G '7/, with T D suppp f 0. 
Because v is piecewise harmonic and [<9 n n]|.B (1+i)i{ nr = *-*’ we 8 et 

ll v MllL2(s (1+ , w \r) = [ S7v-X7(q 2 v) + v 2 \X7q\ 2 dx 

i + ) Jb (1+5)r \t 

= (d n v,q 2 [y 0 v]) + [ v 2 \Vq\ 2 dx. (3.15) 

4S(i+«)i?\r 

We first focus on the surface integral. With the nodal interpolation operator Ih from (13.81) and the 
orthogonality (13.51) . we get 

{d n v,q 2 [q 0 v}) = {d n v,q 2 [q 0 v} - I h (r] 2 hov})) + b(h(v 2 hov]), l) ■ (3.16) 
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The approximation property (13.81 ) leads to 

W^hov] - I h (il 2 [lov))\\ 2 L2{r) < h 2{p+1) Y \ 7 l 2 hov}\ 2 HP+H T) • ( 3 - 17 ) 

TeT h 

Since for each T € Th we have [ 7 o^]|t £ V p , we get D k [-fov]\T = 0 for all multiindices k G No with 
\k\ : = Yli =1 4 = p + 1 and ti\t € "Pi implies 77 1^ = 0 for j € Nq with \j\ > 2 . With the Leibniz 
product rule, a direct calculation (see lFMP13bl Lemma 2] for details) leads to 

\v 2 hov]\ 2 HP+1{T) < ^2 \vh/ov]\ 2 HP{T) + ^4 \hov]\ 2 HP -i {T) , 

where the suppressed constant depends on p. The inverse inequalities (13.101) given in Lemma [331 imply 

Whov] - Ih(V 2 ho v ])\\l HT) < /l 2(p+1) Y (js^\ r l^ v ]\ 2 HP(T) + 

~ jm? II^MII^/ 2 (r) + (^)4 hMWl^B^nnr) ■ (3-18) 

With the trace inequality, we obtain 

Hlov]\\ 2 H i/ 2 (T) = ||7o xt M 

^ IMl£a (n) + II v (^)IIl 2 (q) + IMI L 2 (QC) + IIVMHi 2(nC ) 

- ll' u lli 2 (s ( i + 5 )Jl \r) + ll v (^)lli, 2 (B(i +4 ) Jl \ r ) • ( 3 - 19 ) 

In the same way, the multiplicative trace inequality implies 

h[ 70 ^]|li 2 (r) < ^ Il^lli2 { fl (1+a)li \ r) + ll^llL2 (S(1+5)flXr) hVn|| i2(B(i+i)flVr) . (3.20) 

We apply Lemma [3761 with R = (1 + 5/2)R and 5 = ^75 such that (1 + 5)R = (1 + S)R. Together 
with (13.181) - (13.201) . we get 

\{d n v,r] 2 [^ov} - 4 ( 7 2 [ 7 o^]))| < \\ d nv\\ L 2 {B(i+s/2)Rnr) \\r?[yov] - 4 ( 7 2 M)|| L 2 (r) 

- C (ll Vv lli 2 (4i+7«\ r ) + Jr ll v l^ 2 (4i+7fl\ r )) (ll u llL 2 (i 3 (i + ,) Jl \r) + ll v MII.L2(B ( 1 + 5 )J{ \r) 

+ 44)2 (( 5 r) 1/2 ll u lli 2 (4i+ 5 ) fl\r) + MlL 2 (B( 1 + 5 )fl \r) l^ Vv llL 2 (s ( 1 + 5 )J? \r)) } 

- C JJrJz H Vu lli 2 (4i+*)fl\ r ) + C JJrJi H u lli 2 (S(i + «)fl\ r ) + 4 H v ( r 7 u )llL 2 (s (1+5)fl \r) > 

where, in the last step, we applied Young’s inequality as well as the assumptions ^ < | and 5R < 
2diam(Q) multiple times. The last term in (13.161) can be estimated with (13.81) . 77 < 1, the previous 
estimates (13.181 ) - (13.191) . and the assumption ^ < |, as well as 5R < 2diam(Q) by 

|p(4 (7 2 M),1>| < |p(7 2 [7o^],l)| + \v(v 2 hov] -4(7 2 [7 ou]),1 >| 

^ 41 |4i+a)R nr | 1/2 (|| 7 ? 2 [7 o' u ]|| l 2 ( t ) + II^M — 4(r 2 [7o4)IL 2 (r) ) 

< 41 ((1 + 5)R) {d ~ 1)/2 (||7 2 [7o^]|| L 2 (r) + h 1/2 ||77 2 [7on]||^i /2(r) ) 

< \A ((1 + rn {d ~ m 4lL (B(lw) .\r, + IIvMIIl» ( wot ) ■ 
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Applying Young’s inequality, we obtain 


\v(lh(r) 2 [ 7 ov}).. l)| < C((l + S)R) d 1 \fi \ 2 + IMlia ( B ( 1 + 4 )Jl \r) + j H V Mlli 2 (B ( 1 + 5 )J? \r) ■ 

Inserting the previous estimates in (13.161) . Lemma [3761 Young’s inequality, and the assumption < y 
lead to 

\{dnv,r) 2 [y 0 v])\ < \ (d n v, r] 2 [^ 0 v] - I h {r] 2 [yov})) | + \/i (l h (rj 2 [^ 0 v}), l) | 

- C (YR)2 II Vu ll£ a (fl ( 1 +a)fl \r) + C (^r) 2 ll v lli 2 (S(i + 5 )«\r) + c (( 1 + W * _1 \ e \ 2 

+ 2 H v (^)lli 2 (i3(i +(S ) fl \r) • 

Inserting this in (13.151) and subtracting the term \ |Y(rp;)^ ^ ^, ( ,\iq from both sides finally leads 
to 

ll v ( r ?' (; )llL 2 (B ( 1 + 5 )i? \r) < ^2 ll Vv lli 2 (s ( 1 + 5 ) fl\r) + ^2 ll^lli 2 ( J B(i + 5 )jR \r) + (X 1 + 6 ) R ) d 1 H 2 > 

which finishes the proof. □ 


We consider 7 -shape regular triangulations £n of that conform to Q. More precisely, we will assume 
that every E G £n satisfies either E C 11 or E C O' and that the restrictions £h\vl and £h\q c are 
7 -shape regular, regular triangulations of fl and Q c of mesh size H, respectively. On the piecewise 
regular mesh £h, we define the Scott-Zhang projection .7// : II 1 (W l \ T) —>• Spw := {v : n|o G 
S’ 1,1 |r*) and i;|qc g 5' 1,1 (^//|q c )} in a piecewise fashion by 


J H v = 


Jffv for igU, 
J^v otherwise; 


(3.21) 


here, , J'ff 1 denote the Scott-Zhang projections for the grids £h\q and £h\ci c - Since Jji is a piece- 
wise Scott-Zhang projection the approximation properties proved in HSZ90I1 apply and result in the 
following estimates: 


here, 


V- JHvf H m (E) <CH 2 ^ 


v \ tf*(wg) 


if E c n 

0 < m < i < 1 : 

if E C n c 


(3.22) 


u% = \J{E' e£ H \n : EnE’^Q)}, wf = \J {E' G £h\& : E n E’^ . 


The constant C > 0 in (13.221) depends only on the 7 -shape regularity of the quasiuniform triangulation 
£[[ and the dimension d. 


Let II: ( H l (B R \ T), IHI^) (74, 0 {Br, T p , /r), ||M| \ htR ) be the orthogonal projection, which 
is well-defined since 74, o(Br, T p , / 1 ) C H 1 {Br \ T) is a closed subspace by Lemma [3771 
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Lemma 3.9 Let 5 E (0,1), R E (0,2 diam(f2)) be such that Let B R , B/ 1+s ^ r , B( 1+2 ^ r be 

concentric boxes. Let T p CT be of the form and p E M. Let £ R be an (infinite) 7 -shape regular 
triangulation of M' ; of mesh width H that conforms to Q as described above. Assume -^ < |. Let 
J R : H 1 (R d \ r) —>• 5pw /«' the piecewise Scott-Zhang projection defined in (13.211) . Then, there exists 
a constant C app > 0 that depends only on 12, d,p, and 7 , such that for v E 'Hh,o{B{\+28)Ri T p , p) 

(i) (v — II h,R,[iJHv) | b r € Tihfl{B R , T p , 0); 

(ii) |||u — U hyR}ll J H v\l lhR < C ap p (-^T + -p) ^^3“^ Ill'll/j.,(l+2<5)i? + (U + 2 8)R) ( ' d 1 ^ // “ 

(Hi) dim W < C app , where W := U htR ^ JH'Hh,o(B^ 1+25 ) R , T p , p). 


Proof: For u E Rhfl{B^ l+2 ^ R ,T p , p), we have u E TLh,o{B R ,T p , p) as well and hence 
n h ,R,n (u\b r ) = u\ B r , which gives ©. 

The assumption ^ < | implies [J{E E £h ■ weH B r / 0} C B(\ + j )]R . The locality and the 
approximation properties (13.221) of J R yield 


H 


\ u ^H u llL2( Si? \ r ) + ||V(ti Jhu)Wl2( Br \t) ~ ll^ 7 ' u llL 2 (s (1+5)J? \r) ■ 


We apply Lemma [3781 with R = (1 + 5)R and 6 = y^. Note that (1 + 5)R = (1 + 2 5)R, and < | 
follows from 8 h < 5R = 6II. Hence, we obtain 


^h,R, p JHtr\^ R IIIri/i,it, Jhia) III/!,it — 

= (i?) ~ ^ HU ^L 2 (B R \r) + Jj2 _ ^' u llL 2 (B fl \r) 


< 


H 2 


R 2 


ll Vu llL2(B {1+6)fl \D + 02 ll Vtl ll2 2 (S (1+5 ) ii \r) 


<(*+* 

~ R R 


(1 + 25f 
<5 2 


\\u 


lft.,(i+25)R + (C 1 + 2<5)-R) rf 1 \p\ 2 ^ , 


which concludes the proof dnj). The statement (|TTT]) follows from the fact that 

dim Jii'Hh,o(-B(i+23)it, r p, d) ^ ((1 + 2 5)R/H) d . □ 


Lemma 3.10 Let C app be the constant of Lemma \3 . 91 Let q, k E (0,1), It E (0, 2 diarn(Q)), k E N, 
n/id T p C I he of the form (EH). Assume 


h k q 

R ~ 32k max{C ap p, 1} 


(3.23) 


Then, there exists a finite dimensional subspace W P ofRhjfilf i + k )r, r p , p) with dimension 


dim W k < Cdim 


n + K~ l 

V Q 


d 

k d+1 , 


such that for every v E dlh,o(Bn+ K ) R , F p , p) it holds 
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min ||[7o«] - [ 7 ow]|| L 2 (BHnr ) (3.24) 

w£_Wk 

< C low R(l + n)h~ 1/2 min |u - wf h , 1+lt/2)R 

wEWk 

< c low fi(i + K)h-'i\ k (h», (1+ „ )b + ((1 + K )R )^' 2 Idl) ■ 

The constants C' ( jj m , C\ ow > 0 depends only on il, d, p, and the y-shape regularity of the quasiuniform 
triangulation Th- 


Proof: Let Br and -£>( 1 +sfR with Sj := k( I — ^L) for j = 0,..., k be concentric boxes. We note 

that k = 5q > 5\ >■■■> 5k = y- We choose H = 32kma ^» -jy, where Capp is the constant 

in Lemma [3791 By the choice of H, we have h < H. We apply Lemma [3791 with R = (1 + 5j)R and 
5j = jkfpfjy < Note that 5j-\ = 5j + gives (1 + 5j-i)R = (1 + 2 5j)R. Our choice of H 

implies jL < Hence, for j = 1, Lemma [3~9l provides a subspace W\ of'Hh,o(B(i + s 1 )R>Fp,p) with 
dim W\ < C ^11+ an d a w\ € W\ such that 

l\ v - ' u; l|lk,(l+<5 1 )it ^ 2C'app ^ + ^ J l u IU,(l+5 0 )Jt + (( 1 + do)R) {d 1)/2 

= 8C app — (l + 2Ji) ^IHI^i+^r + ^ (1 + < 5 o) (d 1)/2 

< Q (lHIL,(i+ K )fl + ((! + K )R ) (d ~ 1]/2 \p\) ■ 

Since v — w\ € 'Hh,o{B{i+& 1 )r, r p , 0), we can use Lemma I3T91 again (this time with p = 0) and get an 

approximation W 2 of v — w\ in a subspace W 2 of / Hh,o(B(i+Si)n> r p , 0) with dim W 2 < C ■ 

Arguing as for j = 1, we get 


V-W 1 -W 2 \l h}{ l + 6 2 )R<q\l v 


^l|IU,(l+5i)R ^ V 2 (Mh,(l+K)R + ((! + K ) R ) (d 1)72 h) • 


Continuing this process k — 2 times leads to an approximation w := W/j_| in the space Wj. := 
^2j =1 Wj of dimension dim Wk < Ck ( ^ l+ ^ h ^ = C , Him ((l + n~ 1 )q~ 1 ) d k d+1 such that 


\\v — w\ 


h,(l+K./2)R 


v ~ ™lh,(l+6 k )R < 9* (IN/,,(!+«)* + ((! + K ) R ) (d i)/2 


(3.25) 


The last step of the argument is to use the multiplicative trace inequality. With a suitable cut-off 
function 77 supported by B^ 1+k / 2 )r and 11V 77 11x,°° 1$ (kR)" 1 as well as 77 = 1 on Br, we get for 

2 c H l {B(i +K / 2 )R \ r) 

llbo^ll^^np) < II bo^)]II 2 2 (p) < ll ? ? 2 llL 2 (R d \r)ll 7 7^lli/ 1 ( Rd \ r ) 

~ ^ Z ^LHB ( 1 +k/ 2)r) + 11 Z 11 i 2 (S(i +re /2) ) 11 11 L 2 {B( 1 +K/ 2 ) R \r) 

~ ^ Z ^L 2 (B {1+K/2)R ) + ^ _1 |l 2 lll 2 (S(i +re/2 ) H ) + ^ll Vz ll 2 2(B(i +K /2) fl \r) 

<((1 + k/ 2 )R) h 1 lll 2 lllh,(i+/ C / 2 )J? > 


14 










where the last step follows from the assumption ^ < 1. Using this estimate for z = v — w together 
with (13.251) gives 


min ||[ 70 W] - [low] II L\B R nr) < Ciow(l + n)Rh 1/2 q k + ((1 + n)R) {d 1)/2 M 

w£Wk L 


This concludes the proof. 


□ 


Remark 3.11 The proof of Lemma [3. 101 shows that approximation results in H 1 / 2 can be achieved at 
the expense of an additional factor /i -1 / 2 : With the cut-off function // that is used at the end of the proof 
of Lemma l3.101 we can can bound for z G H 1 (B(i +k / 2 )r \ T) 

H[7oM]||ffi/2 ( r) < IMItfi(R<*\ r) ^ (! + «/2 )Rh~ 1 l4 h , B(1+K/2)R 

Hence, with the spaces IT 7 /,. of Lemma [3.10l one gets 


min ||[7o(r?(n - uj))]|| H i/2 (r) < C[ ow ( 1 + n)Rh 1 q k 

w£Wk 


v IIU,(i+K)i? + ((! + 1)/2 M 


Now we are able to prove the main result of this section. 

Proof of Theorem 13. 11 Choose k = By assumption, we have dist (Br t , Br^) > q 1 diamH/j T = 
\fdq l R r . In particular, this implies 

dist {B( 1+k ) Rt ,B R! t ) > dist (B RT ,B RtT ) - k R T Vd > VdR r ( q" 1 - k) = \fdR r Q - > °- 


Let 4>h G S p,1 (Th ) solve (12.31) . Recall from (12.51) that 


W^hW .ffV 2 (r) + |A| ~ 



l 2 ( r)' 


The potential u = Kfih. then satisfies u G Rh,o(B(i+K)R T i T- A). Furthermore, the boundedness of 
K : LT 1 /2(r) HlJR d ) and -^ < 1 lead to 


Kcf> h 


h,R r ( 1+k) 


< 2 1 + 


R q 


< 


1 + 


Rr 


K<p h 




HRB 2Rt ) 


O’) 


£ ( 1 + 


Rr 


n 12 f 


L 2 (r) 


We are now in position to define the space W^, for which we distinguish two cases. 

Case 1: The condition (13.231) is satisfied with R = R T . With the space provided by Lemma 13. 101 
we set Wk := {[ 70 ^] : w G ITT }. Then, Lemma [3.101 and R r < 2diam(Qj as well as n < 1 lead to 


min 

w&Wk 


— w 


li^B^nr) ~ (1 + n)R T h 1 ^ 2 q k ^ 


Krh 




+ |A| 


< (l+K)(R T + l)h~ 1 / 2 q k n L2 f T2 ^<h~ 1/2 q k n L2 f 

L Z (Y) 


L 2 (r) 
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and the dimension of 144 is bounded by 


dim 144 < C'dim ) kd+1 = Cdim(2 + r\) d q d k d+1 . 

Case 2: The condition (13.2311 is not satisfied with R = R r . Then, we select 444 •’= 
{w | Br nr : w £ S p,1 (Th )} and the minimum in (13.31) is obviously zero. By the choice of k and 
W T > 32k maxfiCapp,!} ’ the dimension of 444 is bounded by 


dim 444 < 



/ 32 A: maxjCapp. 1} 

V 


d -1 




((1 + rj)q- 1 k) d 1 < (2 +r,) d q- d k d+1 . 


This concludes the proof, of the first inequality in (13.31) . The second inequality in (13.31) follows from the 
L 2 (T)-stability of the L 2 (T)-orthogonal projection. □ 


4 %-matrix approximation 

In order to obtain an 72-matrix approximating W 1 (cf. (12.61) 1 we start with the construction of a low- 
rank approximation of an admissible matrix block. 

Theorem 4.1 Fix an admissibility parameter rj > 0 and q £ (0,1). Let the cluster pair (r, a) be 
rj-admissible. Then, for every k £ N, there are matrices X rtr € Mid xr , Y T(T £ MH xr of rank r < 
C'dim (2 + p) d q~ d k d+1 such that 

|| W- 1 ^ - X TtT Y^|| 2 < <7 a p*W (2 ^ 1)/(M ~V- (4.1) 

The constants C apx , C'dim > 0 depend only on 42, d, the y-shape regularity of the quasiuniform trian¬ 
gulation Th, and p. 

Proof: If C dim (2 + p) d q~ d k d+l > min(|r|, |oj), we use the exact matrix block X rcr = W 1 | t xct and 
Y ra = I £ RM X M. 

If Cdi m (2 + rf) d q~ d k d+1 < min(|r|, |oj), we employ the approximation result of Theorem 13. II in the 
following way. Let \ : L 2 (T) —» R be continuous linear functionals on L 2 (T) satisfying A = Sij, 
as well as the stability estimate ||Ai(u;)^i|| L 2 ( r ) < I Mil , 2 (supped for w £ L 2 (T), where the suppressed 
constant depends only on the shape-regularity of the quasiuniform mesh 74 For the existence of such 
functionals, we refer to 1SZ901 . We define ML := {x £ M. N : x t = 0 V i f t) and the mappings 

A r : L 2 (T) —> R T ,r; >-)■ (Aj(u))j 6T and 4> r : M T —>• S p,1 (Th), x Xji/jj. 

j&T 

The interpolation operator 4> r A r is, due to our assumptions on the functionals A,;, stable in L 2 and for 
a piecewise polynomial function f £ S p,1 (Th) we get <J? T (A T ^>) = f |r r with T t := U igT supp C 
Bip. For x £ M r , (12.21) implies 

Ch (d " 1)/2 ||x || 2 < ||$r(x)|| i2(r) < Ch (d " 1)/2 ||x || 2 , Vx € M r . 
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The adjoint AJ : K w —> L 2 (T)' ~ L 2 (T),b i->- Yliex^K °f A z satisfies, because of (12.21) and the 
L 2 -stability of TjAi, 


11 11 z, 2 (r) 


sup 

w£L 2 (T) 


(b, A x w) 2 

IMIz, 2 (r) 


< 


b|| 2 sup 

w&L 2 {T) 


h ( d x )/ 2 ||^ x Azw|| L 2 (r ) 

IMIz, 2 (r) 


< ||b 


Let b € M . Defining / := Ajb| CT , we get bi = (f, ipi) for j 6 <r and supp / C Br^ D T. Theorem l3.ll 
provides a finite dimensional space W k and an element w € Wk that is a good approximation to the 
Galerkin solution fh\B HT nr- It is important to note that the space Wk is constructed independently of 
the function /; it depends only on the cluster pair (r, a). The estimate (12.21) . the approximation result 

from Theorem 13.11 and IT 7 ' 2 / 


< ll/llx 2 (r) < ll A z b lli 2 (r) ~ h {d 1)/2 Il b ll 2 im P!y 


IIA T (f)h ~ A t io| 


< 


< 


l 2 ( r) 

] l -(d-l)/2 _ \ rW 


h -( d -l)/2-l/2 g k 


n L / 


L 2 {V) 


ly<r, < h~WW 
</,-( M -D/2/||b| 


— W\ 


l 2 {b Rt c\T ) 


2 • 


In order to translate this approximation result to the matrix level, let 


W := {A r m : w € W k }. 


Let the columns of X TCT be an orthogonal basis of the space W. Then, the rank of X rcr is bounded by 
dim W k < C'di m (2 + r() d q '“/t 1 . Since X T(T X^ is the orthogonal projection from R ;V onto W, we 
get that z := X Ta XZ a A r d/ ); is the best approximation of A r ^/j in W and arrive at 

\\K<t>h - z \\ 2 < ||A T <f> h - Kw \\ 2 < h-W-W q k ||b || 2 ~ N (2d-m2d-2) q k || b || 2 . (4 . 2) 

Note that A T <f>h = W _1 | T x<r b |cr- If we define Y TiCr := W _1 |^ xo .X rcr , we thus get 2 = X^Y^blo.. 
The bound (14.21 ) expresses 

II (W^lrx, - X rtr Y^) bU|| 2 = ||A r <^ - 2|| 2 < N W-l)/{2d-2) q k || b || 2 . (4.3) 

The space W k depends only on the cluster pair (r, a), and the estimate (14.31) is valid for any b. This 
concludes the proof. □ 


The following lemma gives an estimate for the global spectral norm by the local spectral norms, which 
we will use in combination with Theorem 14. II to derive our main result, Theorem 12.61 


Lemma 4.2 ( IlGraOlll . I Ilac09 , Lemma 6.5.8], HBdrlObl ) Let M € R ;Vx, ' v and P be a partitioning of 
X x X. Then, 


|M || 2 < C s p f y^max{||M | TXff || 2 : (r,cr) € P,level(r) = £} 


\£=0 


where the sparsity constant C S p is defined in ( 12 . 101 ) . 


Now we are able to prove our main result, Theorem 12.61 
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Proof of Theorem 12.61 Theorem 14. 1 1 provides matrices X T(T € M T l xr , Y rcr € Rl <T l xr , so we can define 
the H -matrix by 

w f X r,Y^ if (t, <t) G Pfar, 

w \ W~ 1 | T X(r otherwise. 

On each admissible block (r, <r) G I} :ir we can use the blockwise estimate of Theorem 14. H and get 

IKW 1 - W w )| rXff || 2 < 

On inadmissible blocks, the error is zero by definition. Therefore, Lemma l4~2l leads to 


|W _1 - W-h|| 2 < C sp (^maxllKW- 1 - W w )| rXff || 2 : (r, a) €P,level(r) = £} 


a =o 


< C apx C' sp iV( 2d - 1)/(M - 2) g fc depth(Tx). 

With r = Cdi m (2 + rj) d q~ d k d+1 , the definition b = - l/^j-i) q d ^ d+1 \2 + rj)~ d ^ 1+d ^ > 0 leads 

C A - 


to 


q k = e hvl/(d+1) 5 and hence 


(TV' 1 - W w || 2 < C apx C sp iV( M - 1 )/( M - 2 )depth(T x )e 


_ 6r V(d+i) 


which concludes the proof. 


□ 


5 Stabilized Galerkin discretization 


In the previous section, we studied a saddle point formulation of the hyper-singular integral operator. It 
is possible to reformulate the hyper-singular integral equation as a positive definite system by a rank- 
one correction that does not alter the solution. In numerical computations, this reformulation is often 
preferred, and we therefore study it. Furthermore, it will be the starting point for the "^-matrix Cholesky 
factorization studied in Section [6] below. 

The stabilized Galerkin matrix W st € M ;Vx is obtained from the matrix W € M. NxN as follows: 


w f k = (Wi/j k , ipj) + a (i/j k , 1) (t/jj, 1) = W jk + aB k Bj, Vj, k = 1,..., N. (5.1) 

Here, a > 0 is a fixed stabilization parameter. The matrix W st is symmetric and positive definite. With 
the notation from (12.71) the stabilized matrix W st can be written as 


W st = W + aBEr . 


The interest in the stabilized matrix W st arises from the fact that solving the linear system 


W 


x 


W B 

B t 0 


is equivalent to solving the symmetric positive definite system 


W 


x 


W + aBB T B\ /x 
B t 0 M A 


(5.2) 
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For more details about this stabilization, we refer to f Ste08l Ch. 6.6/12.2]. 

In order to see that the question of approximating (W st ) -1 in the '//-matrix format is closely related to 
approximating W 1 in the "//-matrix format, we partition 


W " 1 




and observe that the inverse (W st ) 1 can be computed explicitly: 

(W st ) _1 = G + (W st ) _1 BP r . 

Hence, the inverse (W st ) 1 can be computed just from a rank one update from G, i.e., a subblock of 
W” 1 . We immediately get the following corollary to Theorem 12.61 


Corollary 5.1 There exists a blockwise rank-(r + 1) approximation 'Wy to (W' ! ) 1 


with 


||(W") _1 - Wft|| 2 < C apx C sp deptHT I )N^ d - i y( 2d - 2 h- brind+1) . 


6 Cholesky decomposition 

In this section we are concerned with proving the existence of a hierarchical Cholesky-decomposition 
of the form W st ~ CyC 1 ^, where Cy is a lower triangular //-matrix. The main results are summa¬ 
rized in Theorem 16.11 It is shown by approximating off-diagonal block of certain Schur complements 
by low-rank matrices. Therefore, the main contribution is done in Section 16.11 the remaining steps 
follow the lines of IIBcb07 1 1GKLB091IFMP13bl . 

The advantage of studying the second system (15.21 ) is that the submatrix W st = W + oBB 7 is sym¬ 
metric and positive definite and therefore has a Cholesky-decomposition, which can be used to derive a 
///-decomposition for the whole matrix. Moreover, the existence of the Cholesky decomposition does 
not depend on the numbering of the degrees of freedom, i.e., for every other numbering of the basis 
functions there is a Cholesky decomposition as well (see, e.g., BHJ131 Cor. 3.5.6]). The existence of 
the Cholesky decomposition implies the invertibility of the matrix W st | pxp for any n < N and index 
set p := {1,... , n} (see, e.g., 1 HJ13 . Cor. 3.5.6]). For the //-Cholesky decomposition of Theorem 16. 1 1 
below we assume that the unknowns are organized in a binary cluster tree Tj. This induces an ordering 
of the unknowns by requiring that the unknowns of one of the sons be numbered first and those of 
the other son later; the precise numbering for the leaves is immaterial for our purposes. This induced 
ordering of the unknowns allows us to speak of block lower triangular matrices, if the block partition 
P is based on the cluster tree Tj. 

The following theorem states that the Cholesky factor C for the stabilized matrix can be approximated 
by a block lower triangular //-matrix and, as a consequence, there exists a hierarchical L //-factorization 
of VV. 

Theorem 6.1 Let W st = CC T be the Cholesky decomposition. Let a partition P ofT xZ be based on 
a cluster tree Tj. Then for every r > 3, there exist block lower triangular, blockwise rank-r matrices 
Cy, L y and a block upper triangular, blockwise rank-r matrix U y such that 

(i) < C' cho iA^d :: idepth(Tx)e~ bri/(d+1) 
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(ii) 


(Hi) 


|W st -C w C w || 

l | w *|| 2 

W-L w u H 


< 2C cho iA rd ^ 1 depth(Ti)e 


_ br i/(d+i) 


+C' chol A fd -idepth(Tx) e 


2 —2br 1 /( d + 1 ) 


w 


< 2(7^01^ d -i depth(Tj)e 


-br !/(d+l) 


+C c 2 hol iV^idepth(T 2 ) 2 e 


2 —2fer 1 /( d + 1 ) 


where C c }, 0 | = C sp C sc y' k, 2 (W Hl ), with the sparsity constant C sp of ( 12. /Oi l, the spectral condition 

, and a constant C sc depending only on il, d, p, the r y-shape 


st-l 


number k 2 ( W st ) := ||W st || 2 W : 
regularity of the quasiuniform triangulation Ty, the admissibility parameter // and the stabilization 
parameter a. 


6.1 Schur complements 


For a cluster pair (r, a) and p := {i € X : i < min(r U <r)}, we define the Schur complement 

S(r, cr) = W st | rXCT — W st | rX p(W st |p X( 0 ) _ 1 W st | pX(T . (6.1) 

As mentioned in HFMP13al such a Schur complement can be approximated by using "H-arithmetic, but 
leads to worse estimates with respect to the rank needed for the approximation than the procedure here. 
Therefore, we revisit our approach from HFMP13all that is based on interpreting Schur complements as 
BEM matrices obtained from certain constrained spaces. 

The main result in this section is Theorem 16.41 below. For its proof, we need a degenerate approxi¬ 
mation of the kernel function k( x,y) = G(x.y) of the single layer operator V given by Vd>(x) : = 
f r G(x, y)(l)(jj)ds y . This classical result, stated here as a degenerate approximation by Chebyshev in¬ 
terpolation, is formulated in the following lemma. A proof can be found in II MP13a l. 


Lemma 6.2 Let rj > 0 and fix f € (0. Iff). Then, for every hyper cube By C W l , d € {2,3} and 
closed D\ C M. d with dist(By , Dy) > rjdia,m(By ) the following is true: For every r £ N there exist 
functions g\y, g 2 y, i = 1,... ,r such that 


K ( x r) -^2,9i,i(x)g 2 ,i{-) 


i— 1 


< c 


( 1 + 1 /^) n i j\~ rl/d 
\d—2 1 ^ " > 


L°°{B y ) 


dist({x}, By) 0 


Vx € Dy, (6.2) 


for a constant C that depends solely on the choice off £ (0 ,2 p). 


The following lemma gives a representation for the Schur complement by interpreting it as a BEM 
matrix from a certain constrained space. A main message of the following lemma is that by slightly 
modifying the Schur complement S(r, a), we can use an orthogonality without the stabilization term. 

Lemma 6.3 (Schur complement and orthogonality) Let (r,a) be an admissible cluster pair, p := 
{i € X : i < min(r U (T )}, and the Schur complement S(r, <r) defined by (16. II) . Let the function 
f € S p,1 (Th ) with (f> = (f> + fip, where f € S p,1 (Th), supp f C T t and fi p € S p,1 (Th ), supp fi p C T p , 
with r T , I f , of the form (13.61) . satisfy the orthogonality 

2 =0 € S p,1 (Th) with suppV> C T p . (6.3) 

Then, there exists a matrix D of rank 2, which is independent off and f, such that 

= (p T (S(r, a) + D)f. 
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Proof: Given 0, 0 is indeed uniquely defined: By definition of 0, we get with the matrix W from (12.71) 


0 = (w0, = (w{0 + 0 P ), ^) i2(r) = (0 T W| rxp + 0jW| pxp )^, 

for 0 € S p,1 (Th), supp '0 C T p and corresponding vector 0 6 Due to /j C /, the matrix W| pxp 
is symmetric and positive definite and therefore invertible. This leads to 

= — 0 T W| rXp W|~x p - 

Thus, we get for 0 with supp 0 <zT a and the vector B from (12.71) that 

(w0,0^+a^ M)GM> = 0 T MTxa + aBB T \ TXa )0 + 0^ {W\ pXrT + aBB T \ pXa )0 

= 0 T (W st | rXCT - W| TXp W|“^ p W st | pX(T ) 0. (6.4) 

With the Sherman-Morrison-Woodbury formula ( IIHJ131 Ch. 0.7.4]), the Schur complement S(r, a) can 
be written as 


S(r,<r) = W st | rX(T - W st | TXp (W st | pxp ) - 1 W st | pX0 - 

= W st | rX(T - (W| rxp + «BB r | TXp ) (W|-^ + P) W st | pX(r , (6.5) 

where P is a rank one matrix given by P = W| P x P «B| p (1 + aB\ T p W\^ p B\ p ) B\ T p W\^ p . Thus, 
comparing the matrices in (16.41 ) and (16.51) . we observe that 

(w0, i/j^+a (fi, {0, 1) = 4> r (S(r, a) + D) 0, 

with a rank-2 matrix D. □ 

Now, we are able to prove the main result of this subsection, an approximation result for the Schur- 
complement S(r, a). 

Theorem 6.4 Let (r, a) be an rj-admissible cluster pair, set p := {i € 1 : i < min(r U <r)}, 

and let the Schur complement S (r, a) be defined in (IQ) . Then for every r > 3, there exists a rank-r 
matrix S r (r, o) such that 


l|S(r, a) - S r (r, <r )|| 2 < C' sc h d ^e~ brmd+1 \ 

where the constants C' sc , b > 0 depend only on Q, d, p, the y-shape regularity of the quasiuniform trian¬ 
gulation Th, and q. Furthermore, there exists a constant C sc depending additionally on the stabilization 
parameter a > 0 such that 

||S(r, a) - S r (r, < t )|| 2 < C sc N ^ d ~^ e -^ 1/(d+1) ||W st || 2 . 

Proof: Let Br t , B /.^ be bounding boxes for the clusters r, a satisfying (12.81) and r p c r defined by 
(13.61) . Lemma l63l provides a representation for the Schur complement as 

0 T (S(r,cr) + B)0 = (w0,0^ ^ +a 1^ (iM)l2(d > ( 6 - 6 ) 
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with the following relation between the functions -?/>, & and the vectors ip, <p, respectively: ip = 
^/=t ''/’jXk 1 where the index j a denotes the j -th basis function coiTesponding to the cluster a, and the 
function <p € S p,1 (Th ) is defined by (p = cp + <p p with (p = ^jXjr an d su PP^p C T p such that 

(w<p,ip\ =0 \hp € S p ’ l (Th) with supp-i/> C T p . (6.7) 

\ / l 2 ( r) 

Our low-rank approximation of the Schur complement matrix S(r, a) will have two ingredients: 
first, based on the the techniques of Section [3] we exploit the orthogonality (16.71) to construct a low¬ 
dimensional space Wk from which for any cp, the corresponding function 6 can be approximated well. 
Second, we exploit that the function ip in (16.61) is supported by T a , and we will use Lemma 16721 

Let 5 = and Bn a , B( ] +( y Rtr be concentric boxes. The symmetry of W leads to 


( W ^) LHT) + “ ($• r> = (a w *) LHn + “ {*• W. 1 W> 

= (*' “'i (W np, ) + (*' W *) + “ (*' *)L 2 (r) ' <6 ' 8) 

First, we treat the first term on the right-hand side of (16.81) . In view of the symmetry property 
S(t, a) = S(< 7 , t) t , we may assume for approximation purposes that diam B Ra < diam B Rt , i.e., 
min{diam(77pj T ), diam(i?R (T )} = \fdR a . Next, the choice of <5 and the admissibility condition (12.81) 
imply 

dist (B (1+2 s)R a ,B RT ) > dist {B R(j 1 B Rt ) - Vd5R a > y/dRa^ 1 - 5) > 0 . 

Therefore, we have (p\B {1+2S)RiT nr p = (7 P | /} n+Mj/Pj nr,, and the orthogonality (16.71) holds on the box 
- 6 ( 1 + 2 5)R a - Thus > b y definition of H h;0 , we have K(p £ 74 j0 (- 6 ( 1 + 25 )^, T P , 0). 

As a consequence, Lemma 13.101 can be applied to the potential Ko with R := (1 + 5)R CJ and k := 
= yyy. Note that (1 + k)(1 + (5) = 1 + 25 and 1 + k -1 = 3 + 77. Hence, we get a low dimensional 

space Wk of dimension dim Wk < Cdi m (3 + rj) d q~ d k d+1 =: r, and the best approximation cp = n ^ (p 
to 6 from the space Wk satisfies 


^ 2 ( s (i+«)no- nr p) 


< R a h~ l > 2 q k 


K<p 


h,{l+25)R CJ 


< ^-l/2 e -bir 1 /( d+1 ) 


tfV 2 (r) ’ 


where we defined b\ := - (3 + v) d ^ l+d ' > > 0 to obtain q k = e &i»- 1/(d+i, _ Therefore, 

C'dim 

we get 


U-^wip) 

< ^-l/2 e -bir 1 /( d+1 ) 

4> 

\ / LRB (1+s)RrT nr p ) 




ffV 2 (r) 


\\WiP\\l^v)- 


(6.9) 


The ellipticity of the hyper-singular integral operator on the screen T p C T, supp(^> — </>) = supp (p p C 
T p , and the orthogonality (16.71) lead to 


< 

m/ 2 (r) ~ 




< 

r%_/ 


||FF(/>|| i y-l/2( r ) 


(p- 


< 

h 1 / 2 (r) ~ 




(p- 


hbr r) 


( 6 . 10 ) 


22 





































Thus, with the triangle inequality, (16.101) . the stability of W : H 1 (T) — y L 2 (T), and the inverse estimate 
(13.91) . we can estimate (16.91) by 




< h~ 1/2 e ~br 1/(d+1) 


L 2 (B ( 1 +s)RtT nr p ) ~ V H^(V) 

< h- 2 e- br ‘ ,U+ " M 11(r| ||i/’ll 1 2 (r) ■ 


+ W&W hV 2 (t) ) 11 11 l 2 (r) 


For the second term in (16.81) . we exploit the asymptotic smoothness of the Green’s function G(-. •). First, 
we mention a standard device in connection with the hyper-singular integral operator, namely, it can be 
represented in terms of the simple-layer operator (see, e.g., I SleOS. Sec. 6]): 




^ curlp <j>, Vcurlp^y , 


( 6 . 11 ) 


where for a scalar function v defined on F, a lifting operator C, and the outer normal vector n, the 
surface curl is defined as 


curlpu = n x 7o nt (V£u), for d = 3, 

curlpu = n • 7o nt (V T £u), V T u = ( 82 V, —d\v) T for d = 2. 


The representation (16.111) is necessary here, since the kernel of the hyper-singular integral operator is 
not asymptotically smooth on non-smooth surfaces F. 

Now, Lemma [6721 can be applied with By = B Ra and D\- = T \ -B(i+< 5 )_r ct , where the choice of 5 
implies 

dist (By,D x ) > * --diam (By). (6.12) 

2y «(1 + Tj) 

Therefore, we get an approximation G r (x, y) = ^' =1 fl\.i( x )fJ 2 .i(y) such that 


||G(x,.)-G r (x,.)|| I » (B ^ ) < Vxer\B (1+J)R „; (6.13) 

here, the constant 62 > 0 depends only on d and r/. As a consequence of (16.121) and (16.131) . the rank-r 
operator W r given by 


<i>, w r i/j 


l 2 (t\b ( 1+s)r<t ) Jt\b (1+s)r 


L 


curlr^(x 


L 


G r (x, y)cmlr'il>(y)ds y ds 1 


% nr 


satisfies with B := (T \ B (x+S)Rci ) x (B Ra D T) 


<t>, (W - W r W 


< 


cur lp 0 


\J meas(r D B Ra 


L 2 (r\B {1+s)RtT ) ~ ‘ ' l 2 (r) 

< h -^s 2 - d R^I 2 e- b2rl/d 


G-Gr 


h V 2 (r) 


IIV’llL 2 (r) 


||curlr^|| L 2 (r) 


< ft-V 1 *’’ 1 " ||0|| t 2 (r) ||^|| l2(r) , 

where the last two inequalities follow from the inverse estimate Lemma [331 the stability estimate (16.101) 
for the mapping 4> >->• (j>, the assumption d < 3 as well as R a < //diam(Q), and the choice 5 = 

Flere, the hidden constant additionally depends on //. 

Since the mapping 


(</>, ij>) Wifi 


+ ($,W r i/>) 

L\B ( 1 +s)RtT nr p ) \ / LRr\B (1+s)RrT ) 


23 





































defines a bounded bilinear form on L 2 (T), there exists a linear operator W r : /- 2 (I j —» L 2 (T) such that 


<j>,Wt!)) + U,W r il)) = (W r (j),ip 

'L\B (1+s)RrT nr,) V’ /i 2 (r\B (1+5)i{(T ) \ 


and the dimension of the range of W r is bounded by 2 r. 
Therefore, we get 




W ^> LH r 


with b := min{6i, 62}. This leads to a matrix S r (r, <r) of rank 2r + 1 such that 

0 T (S(r, cr) + D - S r (r, <r))i/> 


S(r, cr) + D - S r (r, a) 


= sup 
2 0€Kl T l,'i/>GKl CT l 


II0II2M2 


< cfi d - 3 e -^ 1/(d+1) , 


where we have used (12.21) . Consequently we can find a matrix S r (r, a) := S r (r, cr) — D of rank 2r + 3 
such that 

||S (t, a) - S r (r, cr)|| 2 < Ch d - 3 e~ brl/(d+1) . 

The estimate < h~ d+1 (with implied constant depending on a) from I Ste08 , Lemma 12.9] and 

h — N~ 1 ^ d ~ 1 ' ) finish the proof. □ 


6.2 Existence of H-Cholesky decomposition 

In this subsection, we will use the approximation of the Schur complement from the previous section 
to prove the existence of an (approximate) TLCholesky decomposition. We start with a hierarchical 
relation of the Schur complements S(r, r). 


The Schur complements S(r, r) for a block r G Tj can be derived from the Schur complements of its 
sons n, t 2 by 

sfr t) = ^S(ri,ri) S(n,T2) \ 

\S (t 2 ,ti) S(t2,t 2 ) + S(t2,ti)S(ti,ti) 1 S(ri, r 2 ) J 

A proof of this relation can be found in |Beb07| Lemma 3.1]. One should note that the proof does not 
use any properties of the matrix W st other than invertibility and existence of a Cholesky decomposition. 
Moreover, we have by definition of S(r, r) that S(X. T) = W st . 

If r is a leaf, we get the Cholesky decomposition of S(r, r) by the classical Cholesky decomposition, 
which exists since W st has a Cholesky decomposition. If r is not a leaf, we use the hierarchical relation 
of the Schur complements to define a Cholesky decomposition of the Schur complement S(r, r) by 

C(r):= (s(T 2 ,r 1 )(C 1 ( ) r 1 ) T )- 1 C(r 2 )) ’ {6M) 

with S(ti,ti) = C(n)C(ri) T , S(t 2 ,t 2 ) = C(t 2 )C(t 2 ) t and indeed get S(r,r) = C(r)C(r) T . 
Moreover, the uniqueness of the Cholesky decomposition of W st implies that due to CC 7 = W st = 
S (1,1) = C(1)C(1) T , we have C = C(I). 
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The existence of the inverse C(ti 1 follows from the representation (16.141) by induction over the levels, 
since on a leaf the existence is clear and the matrices C(r) are block triangular matrices. Consequently, 
the inverse of S(r, r) exists. 

Moreover, as shown in [ GKLB09 , Lemma 22] in the context of LU -factorizations instead of Cholesky 
decompositions, the restriction of the lower triangular part S(t 2 , ti)(C(ti) t ) -1 of the matrix C(r) to 
a subblock r ' 2 x r{ with t[ a son of r, satisfies 

(S(r 2 ,ri)(C(r 1 ) T )- 1 ) | T / XT / = S(t^,^)(C(^) t ) -1 . (6.15) 

The following lemma shows that the spectral norm of the inverse C(r) - 1 can be bounded by the norm 
of the inverse C(Z) -1 . 

Lemma 6.5 For r € Tj, let C (r) be given by (16.141) . Then, 

max IIC(r) — 1 II„ = ||C(X) _1 || , 

^-rrp || | | Z II I I Z 


Proof: With the block structure of (16.141) . we get the inverse 

r / \-i = ( C(n)- 1 0 \ 

1 \-C(v i )- 1 S{^,T 1 )(C{T 1 ) T )- 1 C{T 1 )- 1 C(r 2 )-V ' 

So, we get by choosing x such that x* = 0 for * € n that 

||C(r ) _1 || 2 = sup [|C(r) _ 1 x || 2 > sup ||C(t 2 ) _ 1 x || 2 = ||C(t 2 ) _1 || 2 

xeRM,|| x || 2 =l " xeKl T 2l,|| x || 2 =l 


The same argument for (C(r) 1 ) 7 leads to 

ll c ( T )“‘ll2 = |( c W _1 ) T || 2 > ll c (n)-‘|| 2 . 

Thus, we have ||C(r ) _1 || 2 > max,; = i :2 ||C(rj ) _1 || 2 and as a consequence max Te Tz ||C'( r ) _1 | 
||C (X) 


-i| 


12 ' 


□ 


We are now in position to prove Theorem 16.11 

Proof of Theorem 16.11 Proof of { 7]): In the following, we show that every admissible subblock t x a 
of C (Z), recursively defined by (16.141) . has a rank-r approximation. Since an admissible block of the 
lower triangular part of C (1) has to be a subblock of a matrix C(r') for some t' € Tj, we get in view 
of (16.151) that C(X)| TXC r = S(r, cr)(C(a) T )~ 1 . Theorem [C4] provides a rank -r approximation S r (r, a ) 
to S (r, a). Therefore, we can estimate 

||C(T)| TXcr —S r (r, cj)(C(cj ) t ) _1 || 2 = ||(S(r, a) — S r (r, a)) (C(cr) T ) _ 1 1 | 2 

< C sc N 2 ^ d -^e~ brl/(d+1) ||(C((T / ) t ) _1 || 2 ||W st || 2 . 

Since S r (r, o){C(o) T )~ l is a rank-r matrix for each ^-admissible cluster pair (r, cr), we immedi¬ 
ately get an Tt-matrix approximation Cy of the Cholesky factor C(X) = C. With Lemma 14.21 and 
Lemma [631 we get 

I1C - C w || 2 < C sc C' sp iV 2 /( d - 1 )depth(T x )e- brl/(d+1) ||C~ 1 || 2 ||W st || 2 , 
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and the equality K 2 (W st ) = K 2 (C) 2 finishes the proof of (0. 

Proof of tfTTTl).- The approximate LIZ-factors Ly, \Jy can be constructed from Gy by 





(6.16) 


where l € solves Cyl = B, and the error estimate follows from 0. 


□ 


7 Numerical Examples 

In this section, we present some numerical examples in dimension d = 3 to illustrate the theoretical 
estimates derived in the previous sections. Further numerical examples about ^-matrix approximation 
of inverse BEM matrices and black-box preconditioning with an decomposition can be found, 

e.g., in I GraOl . Bcb()5b , Gra()5, IB d r 1 Oh 1! FM P13 a 1. where the focus is, however, on the weakly-singular 
integral operator. 

With the choice r/ = 2 for the admissibility parameter in (12.81) . the clustering is done by the standard 
geometric clustering algorithm, i.e., by choosing axis-parallel bounding boxes of minimal volume and 
splitting these bounding boxes in half across the largest face until they are admissible or contain less 
degrees of freedom than ni ea f, which we choose as ni ea f = 50 for our computations. An approximation 
to the inverse Galerkin matrix is computed by using the C++-software package BEM++ liSBA + 15il . 
The %-matrices are assembled using AC A and the C++-library AHMED IIBeb 1211 . 

Our numerical experiments are performed for the Galerkin discretization of the stabilized hyper¬ 
singular integral operator W st as described in Section [5] with a = 1. The geometry is the crankshaft 
generated by NETGEN HSch97H visualized in Figure [Q We employ a fixed triangulation of the 
crankshaft consisting of 5, 393 nodes and 6,992 elements. 



Figure 1: Crankshaft domain 
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Example 7.1 The numerical calculations are performed for the polynomial degree p = 2, resulting 
in N = 13,986 degrees of freedom. The largest block of has a size of 1, 746 2 . In Figure [2j we 
compare the decrease of the upper bound ||l — W st W^|| 2 of the relative error with the increase of 
the block-rank. Figure [3] shows the storage requirement for the computed 71-matrix approximation in 
MB. Storing the dense matrix would need 1,492 MB. We observe exponential convergence in the block 
rank, even with a convergence behavior exp (—hr 1 / 2 ), which is faster than the rate of exp (—hr 1 / 4 ) 
guaranteed by Theorem 12.61 Moreover, we also observe exponential convergence of the error compared 
to the increase of required memory. ■ 




Figure 2: Exponential convergence in block rank 


Figure 3: Exponential convergence in memory required 


Example 7.2 We consider the case p = 3, which leads to N = 31,466 degrees of freedom. The largest 
block of has a size of 3,933 2 . Storing the dense matrix would need 7,608 MB. We observe in 



Figure 4: Exponential convergence in block rank 


Figure 5: Exponential convergence in memory required 


Figure [4] exponential convergence both in the block rank and in the memory. 
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